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1. Introduction 

The limited survey speed of existing radio telescopes means that active galactic nucleus (AGN) 
variability has been studied in detail only for ~ 10^ relatively bright objects, or in only a few epochs 
for larger numbers of sources down to lower flux density limits. However, AGN dominate the radio 
sky at flux densities above a milliJansky (Seymour et al. 2008), and synoptic surveys with SKA 
Phase 1 will provide radio variability information with a cadence of days or less for hundreds of 
thousands of AGN. 

The characteristic timescales of the inner regions of the AGN central engine are accessible 
on human timescales, from hours to days to years, and enable processes occuring close to the su- 
permassive black hole (SMBH) to be probed. Understanding variability on these timescales, and 
correlating variability over a range of radio frequencies with that seen across the electromagnetic 
spectrum, can provide a powerful probe of the connection between disks and jets, inflow and out¬ 
flow, turbulence and shocks, and activity and quiescence. 

In addition to source-intrinsic processes, propagation effects also play an important role in 
observed variability at centimetre wavelengths. Interstellar scintillation (ISS) in the strong refrac¬ 
tive and weak scattering regimes produces apparent intensity variations for sufficiently compact 
sources on timescales of typically days, but sometimes as short as minutes to hours. ISS gives 
access to microarcsecond (ftas) scales (Lovell et al. 2008), allowing the study of changes in intrin¬ 
sic source structure with resolution beyond that achievable with VLBI. Extreme scattering events 
(ESEs) are transient events in which sources undergo large, week- to months-long time-symmetric 
flux density excursions, due to refractive lensing potentially by discrete structures in the interstel¬ 
lar medium (ISM). While these effects have been known for decades, several significant questions 
remain, both with respect to the characteristics of sources that vary and the interstellar structures 
that produce the variation. 

In this Chapter we discuss the contribution that SKA will make to various outstanding ques¬ 
tions in time domain studies of AGN, through observations of both intrinsic and extrinsic variabil¬ 
ity. Although a broad range of specific science outcomes is addressed, similar or identical survey 
strategies are involved, and disentangling intrinsic and extrinsic causes of observed variability will 
be important in order to achieve the desired goals: increased understanding of the physics both of 
AGN and of the intervening media. 


2. Science from time domain studies of AGN: state of the art 


In §2.1 below we discuss direct observations of intrinsic AGN variability, with reference to 


many recent studies. In we discuss variability due to scattering in the interstellar medium, 
including results from the extensive MASIV VEA Survey and follow-up observations, and still- 
mysterious ESEs. Monitoring of ISS probes intrinsic source changes, as well as small-scale struc¬ 


tures in the ISM. In §2.3 we discuss methods for disentangling variations due to intrinsic and 
extrinsic mechanisms. 


2.1 AGN physics from intrinsic variability 

Variability of both total intensity and polarization is a common characteristic of AGN, and 
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occurs with multiple characteristic timescales in individual objects. Variations between different 
bands of the electromagnetic spectrum are often related, with early emission at short wavelengths— 
X-rays or y-rays—eventually propagating to longer wavelengths, including radio, and longer timescales. 
Emission at different wavelengths probes a range of processes occuring in different regions of the 
AGN central engine, and a panchromatic view can help to answer questions of how relativistic AGN 
jets are launched, and how inflow and outflow are coupled. In some cases (e.g. Thyagarajan et al. 
2011) variable radio emission is the only indicator of AGN activity, and can provide a powerful 
probe of buried AGN which show no evidence of an accreting black hole at other wavelengths. 

AGN intrinsic variability exhibits both red and white noise components, as well as quasi- 
periodic oscillations in some cases (King et al. 2013; Wiita 2011), corresponding to variations in 
accretion rate, flares and shocks in disks and jets, transitions between high and low states (KrauB 
et al. 2013), changes in Doppler boosting and jet precession, and other processes taking place 
near the black hole. Long-term monitoring observations are required to sample the full range of 
timescales, to compare to observations at other wavelengths, and to catch rare events such as large 
flares and outbursts. The correlation delay between X-ray and radio can be hundreds of days, so 
long duration campaigns (at least several years) are required for adequate sampling. Day to year 
timescales for AGN variability correspond to the light crossing, viscous, and infall timescales of 
the black hole and accretion disk. Several ongoing programs (e.g. Hughes et al. 1992; Hovatta 
et al. 2007; Lister et al. 2009a; Richards et al. 2011) have monitored a few hundred sources for 
significant periods of time. A few other surveys (e.g. Thyagarajan et al. 2011) have obtained data 
over a handful of epochs for large numbers of sources. But until now, the primary limitation for 
this exercise has been the absence of dedicated radio instruments with the combination of field of 
view and sensitivity to partner with high energy facilities. SKA will open new parameter space 
by providing precision flux measurements as a function of time for almost all of the hundreds of 
thousands of compact sources typical of previous large-area single-epoch surveys. By providing 
high-cadence synoptic monitoring without the potential biases of small sample selections, SKA 
will produce an archive of variability information for a huge range of AGN. 

The launch of the Fermi Space Telescope in 2008 revolutionized gamma-ray studies of AGN 
by providing the first instrument that monitors the entire gamma-ray sky every three hours (Atwood 
et al. 2009). It has become possible to do population studies of gamma-ray loud, and quiet, AGN 
and compare the gamma-ray emission to their radio properties. Lor example, gamma-ray detected 
AGN have larger variability amplitudes in radio (Richards et al. 2011), faster apparent jet speeds 
(Lister et al. 2009b), and higher Doppler beaming factors (Savolainen et al. 2010) than non-detected 
objects. Some studies have also reported that radio and gamma-ray luminosities are correlated 
(Ackermann et al. 2011; Nieppola et al. 2011). All of the above imply a connection between 
radio and gamma-ray variations in the sources. The next generation very high energy Cherenkov 
Telescope Array (CTA), in combination with SKA, will in turn open the possibility to study the 
connection between TeV and radio emission in a large number of sources, probing the acceleration 
of particles up to the highest energies (Giroletti et al. 2014). 

Many existing surveys have focused on blazars, where relativistic beaming and the compact 
sizes of emitting regions produce more variability than in sources oriented further from the line of 
sight. The wide field monitoring capabilities of SKA, however, will enable lightcurves from all 
classes of AGN to be studied, refining unification schemes and our understanding of the geometry 
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of the emitting regions. One caveat is that most, but by no means all (e.g. Kedziora-Chudczer et al. 
2001; Ojha et al. 2010; Stevens et al. 2012), of the existing monitoring surveys were undertaken 
in the northern hemisphere, so not all of the sources from these surveys will be accessible for 
continued monitoring with SKA. 

The details of accretion, and its connection to jet launching in radio-loud sources, as well as 
a wide range of other AGN phenomena, are still not well understood. Radio flares may occur at 
time intervals corresponding to the formation of shocks in the jets, which could be due to interac¬ 
tion with the ISM or torus (Brunthaler et al. 2005), or to changes in accretion rate propagating to 
changes in the jet (Chatterjee et al. 2011). The most commonly used models for flares in blazars are 
so called shock-in-jet models where the variations are caused by shocks propagating down the jet 
and emitting synchrotron radiation (Marscher & Gear 1985; Hughes et al. 1985). Multifrequency 
radio observations can be used to model the shock propagation and jet properties in detail, giving 
insights into the jet geometry and magnetic field (Fromm et al. 2011). Radio flares seen in IIIZw2 
(Fig. ||) and Mrk 348 suggest that jets in Seyferts can accelerate from non-relativistic to relativistic 
speeds during an outburst, but larger samples are required to confirm this (Mundell et al. 2009). 
Interestingly, narrow-line Seyfert I galaxies are an unexpected group of gamma-ray emitting AGN. 
Fermi has detected a handful of NLSl galaxies, confirming that they indeed can also emit high 
energy electrons, and have prominent relativistic jets (Abdo et al. 2009). Modeling these variations 
can probe jet acceleration mechanisms and central engine physics (Kudryavtseva et al. 2011). Un¬ 
derstanding jet production is also critical in studying galaxy formation, since jets play a key role in 
regulating galaxy growth (Hopkins & Elvis 2010). 

Another promising avenue for understanding accretion and jet launching is the study of tidal 
disruption events (TDEs), posited to result from stars that are shredded after passing close to SMBH 
event horizons. The sudden increase in accretion rate results in SMBHs transitioning from a qui¬ 
escent state to a short-lived, typically months-long, AGN-like phase. In some cases radio emission 
is thought to result from the formation of a jet, which interacts with the dense ISM of the galaxy 
(Metzger et al. 2012; Zauderer et al. 2011; Burrows et al. 2011), and produces a radio afterglow that 
peaks around a year after the star is tidally disrupted and lasts for around a year. A large sample of 
TDEs would allow us to address the unanswered question of how jets evolve when accretion rates 
undergo a sudden dramatic increase. By looking for common features in TDEs that are radio loud 
versus radio quiet, their role in AGN radio loudness in general could be constrained (van Velzen 
et al. 2011). TDEs are expected to occur in an average galaxy about once every 10^ yr, so large 
surveys such as those that will be undertaken by SKA are required to see them in significant num¬ 
bers. Eurther discussion of the role of SKA in finding TDEs can be found elsewhere in this volume 
(Donnarumma et al. 2014). 

A recent study (MacEeod et al. 2010) found a dependence of AGN variability on black hole 
mass, luminosity, and Eddington ratio, as well as larger variability amplitudes for radio-loud ob¬ 
jects than for radio-quiet. However, the relationship between variability and the properties of the 
host galaxy and black hole, and the accretion dynamics, are by no means clear (Zuo et al. 2012). 
Recently, it was shown that the jet magnetic field and accretion disk luminosity are tighly corre¬ 
lated in blazars (Zamaninasab et al. 2014). This result was obtained by comparing black hole mass 
estimates to the jet magnetic flux, calculated using radio interferometric observations. This result 
shows that it is possible to obtain information about the detailed physics near the black hole, by 
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Figure 1: Lightcurve of III Zw 2 at 22 and 37 GHz from the Metsahovi Radio Observatory, showing flaring 
activity suggestive of acceleration of jets from non-relativistic to relativistic speeds during outbursts. Credit; 
Terasranta et al., A&A, 440, 409, 2005, reproduced with permission © ESO. 


relating observations further out in the jet to the smallest scales near the black hole. 

VLBI observations provide direct evidence for relativistic outflows from AGN, and have been 
used to monitor the motion of jet components on parsec scales (e.g. Lister et al. 2013). It is possible 
to determine jet compositions, magnetic fields, variations in jet orientation, ages and densities of 
components, interplay between the jet and ambient medium, as well as other physical conditions 
in the jet and close to the black hole by comparing parsec-scale variations with spectral energy 
distributions, polarization and Faraday rotation, as well as observations at other wavelengths and 
relativistic MHD simulations (Hovatta et al. 2012, 2014). VLBI components are by their very 
nature intrinsically variable (Lister et al. 2009a), but are also subject to variability due to interstellar 
scattering (§|2.2|). Flares and superluminal blobs may be connected to changes in magnetic field 
on timescales of years, or to current-driven instabilities, tangled magnetic fields, or reconnection 
events, and SKA can help probe this (Agudo et al. 2014). 

AGN activity appears to be enhanced during galaxy mergers (Manzer & De Robertis 2014) 
and in overdense environments (Croft et al. 2007). However, the relationship between mergers and 
the triggering of various kinds of AGN, including those selected at different frequencies, is still 
not clear (Kocevski et al. 2012; Karouzos et al. 2014). The details of what happens to SMBHs 
after mergers are also not yet well understood. In some instances, AGN pairs (or multiple AGNs) 
with an offset in position and / or velocity are seen in galaxies at various stages of mergers (e.g. 
Komossa et al. 2003; Comerford et al. 2011). With large area, sensitive, high resolution SKA 
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Figure 2: Centimetre-band lightcurves of BL Lac object OJ 287 from the University of Michigan Radio 
Telescope. Optical lightcurves of this source, including archival data dating back over 100 years, show a 12 
year periodicity interpreted as the orbital period of a secondary SMBH intercepting the accretion disk of the 
primary. The radio lightcurves shown here (from Aller et al. 2014, reproduced by permission of the AAS) 
show flaring activity on shorter timescales. From top to bottom, the panels show total flux density, fractional 
linear polarization, and polarization electric vector position angle. 


surveys, many more dual AGN (including binary AGN, where the AGN are clearly orbiting each 
other, e.g. Rodriguez et al. 2006; Deane et al. 2014) or systems containing an AGN and a quiescent 
SMBH (e.g. Comerford & Greene 2014), should be detected at various stages in their evolution. 
In at least one case known to date, the blazar OJ 287, variability of optical emission is proposed 
to arise as a secondary black hole plunges through the accretion disk of the primary (Valtonen 
et al. 2008). The optical lightcurve shows periodic outbursts every 12 years, corresponding to the 
proposed orbital period, and radio monitoring (Fig. ^ shows evidence of long term cycles of radio 
activity as well as periodic variability associated with flares and jet features (Hughes et al. 1998). 

For closer binaries, the “final parsec problem”—the ineffectiveness of dynamical friction with 
stars near SMBHs in shrinking their separation to < 1 pc (Milosavljevic & Merritt 2001)—remains 
unsolved, and likely requires a more detailed understanding of accretion disk dynamics. Periodic 
variability due to jet precession, a signature of gravitational wave dominated inspiral of SMBHs 
during the final stage of mergers (Haiman et al. 2009), can be used to search for merging systems 
containing a radio-loud AGN. A detection of this signal, which should be within reach of SKA and 
its precursors (O’Shaughnessy et al. 2011), could be used to probe the physics of viscous orbital 
decay as a function of black hole mass, and confirm that inspiral is driven by gravitational wave 
emission. 

By probing the inner regions of AGN central engines, studies of AGN intrinsic variability 
with SKA will inform two Key Science Projects: “Strong Field Tests of Gravity Using Pulsars and 
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Black Holes”, and through the connections between infall and outflow in AGN, and feedback on 
surrounding gas, “Galaxy Evolution, Cosmology, and Dark Energy”. 

2.2 Science from interstellar scattering 

Radio wave propagation effects and implications for SKA, applying to studies of both pul¬ 
sars and AGN, were discussed in some detail by Eazio et al. (2004). Here we provide an updated 
overview specifically related to AGN variability due to interstellar scattering, and outline the cur¬ 
rently outstanding questions that can be addressed with SKA. 

2.2.1 Interstellar scintillation 

Over the past two decades, evidence has accumulated to demonstrate that rapid—intra-day and 
inter-day—flux density variations of AGN at centimetre wavelengths are predominantly caused 
by ISS rather than intrinsic variability. This evidence includes the frequency dependence of the 
observed variability (Kedziora-Chudczer et al. 1997), measurement of variability pattern arrival 
time delays of several minutes between widely separated telescopes for the most rapid variations 
(Jauncey et al. 2000; Dennett-Thorpe & de Bruyn 2002; Bignall et al. 2006), annual cycles in 
the characteristic timescale of variability (e.g. Dennett-Thorpe & de Bruyn 2001; Rickett et al. 
2001; Jauncey & Macquart 2001; Qian & Zhang 2001; Bignall et al. 2003) and a strong correlation 
between cm-wavelength intraday variability (IDV) and Galactic Ha emission measure along the 
line-of-sight, showing that this IDV is related to the ionized Galactic ISM (Eovell et al. 2008; Koay 
etal. 2011). 

Observed ISS implies the presence of high brightness temperature compact components, no 
larger than tens of /ras in angular size (e.g. Eovell et al. 2008). ISS, with an adequate model of 
the scattering parameters, can thus be used to estimate source angular size and compact fraction, 
as well as investigate frequency- and polarization-dependent sub-structure, at resolutions beyond 
those achievable with ground-based VEBI. Angular scales of this order are presently being probed 
directly on baselines to the RadioAstron Space Radio Telescope (Kardashev et al. 2013), but the 
sensitivity is limited to bright sources. Studies of ISS of fainter sources over wide bandwidths 
probe a fundamentally different, much weaker, regime of AGN jet power. 

The MASIV VEA Survey at 5 GHz (Eovell et al. 2003, 2008) found that more than half of 
all compact, flat-spectrum sources > 100 mJy exhibited rms variations > 2% on timescales < 3 
days due to ISS, with increased fractional variability towards lower flux densities. This trend is as 
expected for brightness temperature-limited sources, where angular size is expected to scale as s/S 
assuming no change in the Doppler boosting factor. The majority of MASIV sources have firm 
spectroscopic identifications—all AGN, with almost 80% being flat-spectrum radio quasars and 
most of the remaining fraction BE Eac objects (Pursimo et al. 2013). 

A suppression of ISS towards higher redshift sources (Eovell et al. 2008) has been proposed as 
a possible indicator of angular broadening due to scattering in the turbulent intergalactic medium 
(IGM), making the high redshift sources appear too large to scintillate through the Galactic ISM. 
However, a careful analysis of the MASIV Survey and follow-up VEA data at 5 and 8 GHz found 
no conclusive evidence for IGM scatter-broadening of the high redshift MASIV sources. After 
accounting for a steepening with redshift of average spectral index between 5 and 8 GHz, the 
decrease in amplitude of ISS with redshift was found to be no longer significantly in excess of 
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that expected if angular size 6 oc VT+z, as expected for a brightness temperature and flux limited 
sample of sources (Koay et al. 2012). 

The statistics of short timescale radio variability for a larger and much fainter source sample, 
obtainable with SKAl, will allow tighter limits to be placed on extragalactic scatter-broadening 
through the suppression of ISS. As lower flux density AGN cores can have smaller intrinsic source 
sizes, they are more likely to be dominated by scatter-broadening. Observations of large num¬ 
bers of faint sources in fields at high Galactic latitudes would provide a means of determining the 
ISM contribution to scattering; source redshifts are also required in order to measure, or obtain 
limits on, intergalactic scattering. Koay & Macquart (2014) have recently shown that while scatter¬ 
broadening in the warm-hot component of the IGM is not expected to be detectable, intergalactic 
scatter broadening can be significant, ~ 100/ras at 1 GHz and ~ 3/ras at 5 GHz, for sight-lines in¬ 
tersecting within a virial radius of at least one galaxy halo. Targeted, multi-frequency ISS surveys 
of compact sources behind rich clusters, galaxies or Damped Lyman-Alpha Systems may provide 
the best chance of detecting extragalactic scatter-broadening through the predicted angular size 
scaling as providing a unique probe of the turbulent properties of these regions. Evidence 
for scatter-broadening in an intervening galaxy has been found from VLBI observations of AGN 
behind M31 (Morgan et al. 2013). 

A handful of sources have been found which exhibit unusually rapid variability on timescales 
< 1 hour (Kedziora-Chudczer et al. 1997; Dennett-Thorpe & de Bruyn 2000; Bignall et al. 2003), 
which is modelled as ISS in the weak scattering regime due to very nearby scattering plasma, 
within tens of parsec of the Sun (Rickett et al. 2002; Dennett-Thorpe & de Bruyn 2003; Bignall 
et al. 2006). The scattering is determined to be highly anistropic, as has been found also for a num¬ 
ber of pulsars displaying diffractive “scintillation arcs” in their secondary spectra (Walker et al. 
2008; Brisken et al. 2010, and references therein). This suggests scattering in thin “screens” in 
which the electron density fluctuations are subject to extreme magnetic stresses (Tuntsov et al. 
2013). The specific nature of the material responsible is still unknown. Despite two of the three ex¬ 
emplary “intrahour variables” being discovered serendipitously, the MASIV Survey found no more 
examples of such rapid, large-amplitude ISS, implying that the covering factor of nearby screens is 
very small; <C 1% of compact sources would be expected to vary on intra-hour timescales. 

The correlation between H-alpha and short-term variability from MASIV suggests that in gen¬ 
eral, the more typical slower, low-amplitude ISS is probably related to the overall amount of ionized 
material in the ISM along the line-of-sight, but this is not the case for the intrahour variables. Be¬ 
cause AGN angular sizes are typically larger than the Fresnel scale at interstellar distances, they 
are “resolved” by ISS. The smaller the screen distance L from the observer, the larger the angular 
Fresnel scale, so nearby scattering will tend to dominate the observed variability; not only is the 
ISS less quenched, but nearby screens also induce faster variations compared with more distant 
scattering material, since the Fresnel zone crossing time which sets the timescale for weak ISS 
scales as s/L for a given scintillation velocity. ISS behaviour is often observed to be episodic, 
which can be due to a combination of intermittency in interstellar turbulence and intrinsic changes 
in source structure, although the relative importance of the two effects remains poorly constrained. 
This problem can be addressed by comparing longer term source-intrinsic changes in spectral in¬ 
dex, flux density, polarization and VFBI structure, with changes in ISS behaviour for a large source 
sample. 
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Figure 3: The extreme scattering event toward QSO 0954+658, the prototypical and exemplar ESE. The 
top panel shows the 2 GHz flux densities measured by the Green Bank Interferometer, and the bottom panel 
shows the 8 GHz flux densities. Data are from Rickett et al. (2006). 


Finding a large number of rapid intraday variables through SKA Surveys will help to determine 
the nature of the local scattering plasma in the ISM. While fast intra-day variability can be picked 
up in undersampled data, extracting useful science would require detailed follow-up observations, 
either with (a sub-array of) SKA or with other instruments. Well-sampled ISS light-curves across 
a broad bandwidth, ideally several GHz, can be used to determine the detailed polarization- and 
frequency-dependent structure of the ultracompact jet with /ras precision. This information can be 
used to model the magnetic field configuration, geometry, particle density and pressure in the jet 
(Rickett et al. 2002; Macquart et al. 2013). 

2.2.2 Extreme Scattering Events 

Extreme scattering events (ESE) are a class of dramatic changes in the flux density of radio 
sources (Eiedler et al. 1994). They are typically marked by a decrease, as large as > 50%, in the 
flux density near 1 GHz for a period of several weeks to months, bracketed by substantial increases, 
viz. Eigure Because of the simultaneity of the events at different wavelengths, the time scales 
of the events, and light travel time arguments, ESEs are likely due to strong scattering by ionized 
structures in the Galactic ISM (Eiedler et al. 1987a; Romani et al. 1987, but see Walker & Wardle 
1998). Eirst identified in the light curves of extragalactic sources, ESEs have since been observed 
during a timing program of the pulsars PSR B1937+21 (Cognard et al. 1993; Eestrade et al. 1998) 
and PSR 11643-1224 (Maitia et al. 2003). 

While the focus of this chapter is on the variability of AGN, ESEs probably are related to 
“strong fringing events” observed in dynamic spectra from pulsars (Hill et al. 2005). In those 
events, the dynamic spectra—intensity as a function of time and frequency—show “fringes,” con¬ 
sistent with the interference from multiple images of a pulsar resulting from strong lensing along 
the line of sight to the pulsar. 

Modeling of ESE light curves leads to inferred densities rig >10^ cm^^ within these ionized 
structures (Romani et al. 1987; Clegg et al. 1998). In turn, these densities imply pressures nT ~ 
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10^ K cm ^ or more, well in excess of the “average” interstellar pressure nT ~ 3000 K cm ^ 
(Kulkarni & Heiles 1988). 

Outstanding issues related to ESEs include the following: 

• What is the mechanism by which they occur? Most models predict that ESEs arise from 
strong refraction through relatively isolated structures. Observational support for this pre¬ 
diction is mixed. Pulse times of arrival from the pulsars that have undergone an ESE do 
show an increased scatter, but no increase in the pulse width—as would be expected from 
the extra path length caused by the refraction. However, VEBI observations of the source 
B1741—038 while it was undergoing an ESE show its angular diameter to increase, contrary 
to expectations from a purely refractive model (Eazio et al. 2000). One obvious resolution 
is to model the ESE structures as being internally turbulent, but largely refracting struc¬ 
tures. The serendipitous detection of multiple imaging of the quasar 2023-1-335 observed 
with VEBI during an ESE was recently reported, and found to be consistent with strong 
refraction through a turbulent plasma structure (Pushkarev et al. 2013). 

• What is the relationship of ESE structures to other phases of the ISM? or are ESE struc¬ 
tures even interstellar? ESE structures may represent relatively isolated structures, perhaps 
in pressure balance with a lower density, higher temperature “background” phase (Clegg 
et al. 1988); they may represent manifestations of the most extreme pressure ranges within 
the ISM and thereby trace energy input into the ISM (Jenkins & Tripp 2001); they may re¬ 
sult from a low level of cosmic-ray ionization within an otherwise neutral tiny scale atomic 
structure (TSAS; Heiles 1997); or they may not be interstellar at all, but due to photoionized 
molecular clouds in the Galactic halo (Walker & Wardle 1998). HI absorption measurements 
during the ESE of B1741—038 show no statistically significant change in the HI line during 
the ESE (Eazio et al. 2001a). However, the limits on the changes in optical depth, and there¬ 
fore HI column density, only marginally rule out a connection between ESE structures and 
TSAS in the ISM. Any change in the HI column density is limited to be ANh < 10^^ cm^^, 
while TSAS has a typical column density of 3 x 10^^ cm^^. Geometrical factors (e.g., the 
line of sight did not pass through the center of the cloud) could easily account for this differ¬ 
ence. Conversely, the velocity range of the HI absorption measurements was no more than 
250 km s^^ If ESE structures are halo objects, the allowed velocity range could approach 
1000 km s^^ Thus, the existing Hi absorption observations cannot constrain significantly 
either of the two competing models. 

• What is the frequency with which ESEs occur and what is their sky distribution? The Green 
Bank Interferometer (GBI) monitored approximately 30 sources for over a decade and an 
additional 120 sources for five years. A fofal of 12 ESEs were found in fhis moniforing 
program. However, a six-monfh moniforing program wifh fhe NRAO 300 ff felescope found 
a comparable number of ESEs (R. Eiedler 1998, privafe communicafion). The 12 ESEs 
observed in fhe GBI program appear fo be locafed near fhe edges of radio Eoops I-IV (Eiedler 
el al. 1994), sfrucfures fhoughf fo be old supernova remnanfs. If Ibis is correcf, old supernova 
remnanfs could provide fhe high-pressure environmenf in which ESE lenses could survive 
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(Clegg et al. 1988). However, the radio Loops cover a large fraction of the sky, and the sky 
distribution might also be consistent with a halo population of objects. 

Since the cessation of the GBI monitoring program in 1994 (Fiedler et al. 1987b; Waltman 
et al. 1991; Lazio et al. 2001b), there has been no monitoring program observing a sufficiently 
large number of sources, at a low enough radio frequency, and with frequent enough time sampling 
to have a realistic chance of detecting the requisite number of ESEs to address these questions. 
Although some new programs aim to tackle this in the lead-up to SKA Phase 1, e.g. the ASKAP 
VAST project (Murphy et al. 2013), these cannot match the efficiency of potential SKA surveys to 
detect a sufficient number of ESEs to address their statistics (§^4]). 

Identification of ESEs as they occur is critical in order to determine the nature of the refracting 
“lenses” responsible. Multiple imaging on milliarcsecond scales is predicted at strong caustics, 
and VEBI monitoring of ESEs will differentiate between competing models for ESEs (Pen & King 
2012). Detailed radio follow-up over a broad frequency range, and observations at other wave¬ 
lengths, including optical/infrared, will also help to rule out or support various models. 

2.3 Determining the origin of AGN radio variability 

Observations at multiple frequencies and classification of light curves based on their structure 
function, power spectral density, and other statistics will help to disentangle intrinsic and extrinsic 
causes of variability. Ofek & Erail (2011) argued that a large fraction of variability at 1.4 GHz 
must be RISS, based on a flat structure function between ~ 200 — 1800 days for a sample of 43 
variables. Rickett et al. (2006), in an analysis of GBI monitoring data for 146 sources, applied a 
simple filter on the time series to separate slow intrinsic variations from fast ISS, confirming that 
both processes are important at the observed frequencies near 2 and 8 GHz (see Eigure 

It is expected that propagation effects may dominate at low frequencies, below ~ 2 GHz, while 
large-amplitude variations at higher frequencies, > 8 GHz are more likely intrinsic — at least on 
timescales of weeks to months and longer, although both ISS (Savolainen & Kovalev 2008; Kara 
et al. 2012) and ESEs (Pushkarev et al. 2013) have been observed at 15 GHz. However, a test 
shows no statistical correlation between line-of-sight Galactic Ha emission and 15 GHz variability 
on timescales of months to years from published OVRO data, indicating intrinsic variability domi¬ 
nates, whereas there is a strong correlation between Ha emission measure and inter-day variability 
at 5-8 GHz, indicating that the latter is predominantly ISS. At these intermediate frequencies, both 
intrinsic and ISS-induced variability are likely to be important. Concurrent lightcurves, includ¬ 
ing polarization information, for large numbers of AGN, using Band 5, 4.6-13.8 GHz, of SKAl- 
MID—and MeerKAT’s X-band receivers—along with two lower frequency bands, will provide a 
powerful discriminant between intrinsic and extrinsic variability. AGN lightcurves from SKAl- 
EOW, at 50-300 MHz, and SKAl-SUR, at 650-1670 MHz (PAP Band 2), will likely be dominated 
by extrinsic variability, although the panchromatic view provided by the wide frequency coverage 
of SKA, in conjunction with instruments such as the Earge Synoptic Survey Telescope (ESST) and 
high energy facilities, will provide a powerful discriminant of the source of variability for a range 
of objects. 

The largest amplitude variations due to ISS are usually seen close to the transition between 
weak and strong scattering, which is dependent on Galactic latitude but typically ~ 5 GHz (Kedziora- 
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Figure 4: Flux density at 2 and 8 GHz for 1514+197 from the Green Bank Interferometer monitoring 
program, analysed by Rickett et al. (2006) (figure reproduced with permission). The solid curve is the time 
series (Jy) smoothed over 365 days with a truncated cosine function. It is interpreted to be dominated by 
intrinsic variation in the synchrotron radiation from the source. The short-term deviations from the line are 
due to ISS. The trace below each curve gives the rms error in a single measurement (on the same scale) 
relative to the horizontal axis. This demonstrates the value of long-term, multi-frequency monitoring with 
high cadence to determine the origin of variability. 


Chudczer et al. 1997; Walker 1998). When the source is much larger than the angular size of the 
first Fresnel zone at the scattering screen distance, then weak ISS will be suppressed, but slower 
refractive ISS may still be observed at lower frequencies in the strong scattering regime. Intrinsic 
variability at centimetre wavelengths will typically have longer timescales and larger amplitude 
than scintillation. A preliminary comparison of modulation indices from four years of OVRO data 
(Richards et al. 2014) at 15 GHz — predominantly intrinsic variability — and short-timescale (2 
day) modulation indices from the MASIV Survey at 5 GHz — predominantly ISS — for sources 
common to both samples shows a significant correlation: sources which show the largest amplitude 
ISS are also more intrinsically variable. This finding, yet to be published, indicates that intrinsic 
source compactness, core dominance and beaming are important factors for ISS, in addition to the 
properties of the ISM on a particular line of sight. 

ESEs may be distinguished from intrinsic variability or “normal” ISS through their strong 
frequency dependence, the presence of caustics at frequencies of at least a few GHz, and the overall 
time-symmetric behaviour in the light curves, including a sharp flux density decrease (see Eigure|^. 
As discussed in Section 2.2.2| , identifying ESEs as they occur is critical to address the nature of 
the leasing structures responsible, by obtaining detailed follow-up observations at high angular 
resolution and at other wavelengths. 

Observations of variability, both intrinsic and due to ISS, may additionally provide a means 
of distinguishing between emission due to AGN and star formation for large samples of the faint 
radio source population. The survey speed of SKA will enable high cadence monitoring of sample 
sizes vastly in excess of what is achievable with existing telescopes. 
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3. Time domain studies of AGN with the SKA 

3.1 Source statistics - expected rates of variability 

The VLA FIRST Survey (Beeker et al. 1995) measured an average souree density of ~ 90 per 
square degree down to 1 mJy at 1.4 GHz, the vast majority of whieh are AGN. Below this flux 
density, star forming galaxies heeome more numerous, although a signifieant fraetion of sourees 
are found to be AGN-driven (see e.g. Figure 4 of Norris et al. 2011, and referenees therein). The 
mJlVE-20 VLBA Survey of FIRST sourees (Deller & Middelberg 2014) found an inereasing frae¬ 
tion of eompaet eomponents down to 1 mJy, with 30-35% of sourees being eore-dominated and 
~ 10% being eompletely eompaet. Sueh eompaet sourees are likely to be variable. Deller & Mid¬ 
delberg (2014) obtained a lower limit on the highly variable souree fraetion at 1.4 GHz based on the 
number of sourees whieh showed an inerease in the VLBl flux density over the VLA FIRST Sur¬ 
vey flux density measured some 10-15 years earlier. This number is an underestimate of the true 
number of variables beeause the different resolutions of the VLA and the VLBA mean that only 
a relatively large inerease in flux density of eore-dominated sourees eould be measured, sinee an 
observed deerease in flux density eould be entirely due to resolution effeets rather than variability. 
Based on this, a eonservative lower limit to the variable fraetion at 1.4 GHz, for > 25% variability 
on timeseales of years, is 0.7% of all radio sourees in the 1-20 mJy range. This is in agreement 
with results from eomparing flux densities for sourees deteeted in multiple epoehs of the LIRST 
survey, where around 0.6% of sourees vary by 25% or more (Thyagarajan et al. 2011). Down to 
flux density levels of 40 /tJy, Mooley et al. (2013) found 1% of unresolved sourees variable at the 
4a level on timeseales from 1 day to 3 months, with no evidenee that the fraetional variability 
ehanges along with the known transition of radio-souree populations below 1 mJy. Optieal identi- 
fieations showed the variable radio emission to be assoeiated with the eentral regions of an AGN 
or star-forming galaxy. 

The fraetion of sourees identified as variable naturally inereases as souree lighteurves are mea¬ 
sured with higher preeision and over a wider range of timeseales and frequeneies. Close to the 
deteetion threshold of any survey, it is diffieult to determine whether a souree is variable or not 
(see, e.g. fig. 26 of Croft et al. 2013, for a quantitative example for the ease of the Allen Teleseope 
Array). A strawman survey with SKAl-SUR or SKAI-MID in band 2, ~ 1.2 GHz, would be eapa- 
ble of surveying the entire visible sky every day to an RMS sensitivity of around 100 /tJy beam^^, 
enabling 10a detentions of variability at the 50% level for sourees brighter than ~ 2mJy. Almost 
all of these sourees will be AGN. Of order 10^ sourees, similar to the entire eatalogs of previous 
surveys sueh as LIRST and NVSS, eould be monitored for sueh variability on a daily basis, or 
larger numbers of fainter sourees by binning data. Lor bright sourees, high-preeision lighteurves in 
total as well as polarized flux density will enable measurements of variability with an order of mag¬ 
nitude or more better preeision than is possible with existing faeilities. The aeeuraey and stability 
of ealibrated flux density measurements is a eritieal parameter here; as an example, the MASIV 
VLA Survey aehieved <1% rms stability over eaeh observing session. Beeause of the large num¬ 
ber of available in-beam ealibrators, many of whieh will be steep-speetrum and not signifieantly 
variable, in prineiple SKA instruments should be able to aehieve at least this level of preeision for 
flux density measurements, but eareful implementation of ealibration proeedures will be required. 
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At flux densities < 1 mJy, the population begins to transition from AGN-dominated to starburst- 
dominated, although even at 1.4 GHz flux densities of ~ 50/iJy, AGN still make up ~ 25% of radio 
sources (Seymour et al. 2008). Starburst galaxies are not expected to be variable at radio wave¬ 
lengths, with the exception of rare transient events such as radio supemovae, but low-luminosity 
AGN—including those that are technically classified as “radio-quiet” but nevertheless have low- 
power radio emission—will still be present, both in ordinary galaxies and in some of the starburst 
galaxies themselves. The high resolution—tens of milliarcseconds—enabled by SKA2 will help 
to distinguish nuclear point sources from emission due to a surrounding disk, but variability will 
also be an important clue to the presence of AGN whose emission might otherwise be swamped by 
the host galaxy. Variability may also be key in identifying and classifying counterparts to sources 
found by CTA or other high energy facilities, where positional coincidence alone may be insuffi¬ 
cient (Giroletti et al. 2014). 


3.2 Frequency coverage 


Intrinsic AGN variability is more extreme towards higher frequencies. Therefore SKAl-MID 
band 5, covering 4.6-13.8 GHz, along with two lower frequency bands, would be extremely useful 
for distinguishing intrinsic from extrinsic changes. While not critical for this work, an expansion 
to 24 GHz as part of SKA2 would enable monitoring observations at similar frequencies to existing 
long-term surveys focusing on intrinsic variability (e.g.. Fig. [I]). 

The inclusion of SKAl-MID or SKAl-SURband 2, covering 1.4 GHz, would allow measure¬ 


ments of HI absorption towards sources undergoing ESEs (section [2.2.2D . The optimal frequency 
range for ESE searches is not yet well determined. Refractive effects are stronger towards lower 
frequencies, and the “cross-section” for events occurring may also be higher due to larger source 
sizes, however larger source sizes are also a problem as they result in the ESE light curves being 
“smeared out”. ESEs have also been seen at frequencies as high as 15 GHz (Pushkarev et al. 2013). 
Ideally, monitoring would be done over a wide range of frequencies; this could be facilitated by 
the use of sub-arrays observing simultaneously in different frequency bands. The appearance of 
caustics at frequencies of a few GHz helps to unambiguously identify ESEs (see Eigure IH). 

Eor ISS, the largest amplitude variations may be observed close to the transition frequency 
between weak and strong scattering, typically a few GHz. Weak ISS occurs on shorter timescales 
than refractive ISS. In general, AGN components are too large to exhibit the narrow-band, short- 
timescale diffractive ISS observed for pulsars in the strong scattering regime (Dennison & Condon 
1981, but see Macquart & de Bruyn 2006). 


3.3 Modeling intrinsic variations 

In the brightest AGN, blazars, the variability is often modeled with shock-in-jet models (Marscher 
& Gear 1985; Hughes et al. 1985). Recently, detailed modeling of blazar variability during high ac¬ 
tivity periods was conducted using 4.8, 8 and 15 GHz total intensity and linear polarization archival 
data from the University of Michigan Radio Astronomy Observatory (Aller et al. 2014). The vari¬ 
ations were modeled with shocks moving down the jet, and the results showed that no extreme 
jet conditions are required to form the shocks. The jet and shock parameters obtained were also 
very similar in all sources modeled. Multiple frequencies providing spectral index information and 
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inclusion of linear polarization were essential for constraining the model parameters. The higher 
frequencies of SKAl-MID will provide information about the intrinsic variability of a much larger 
set of AGN, and modeling of a large sample of sources will become possible. In particular, SKA 
will have unprecedented ability to monitor changes in the weak, typical V jl < 1%, but highly vari¬ 
able circular polarization (CP) component (Homan et al. 2006), a powerful probe of the evolution 
of magnetic fields in relativistic jets (Agudo et al. 2014). CP is associated with the most compact 
components of AGN, which are also subject to ISS. Unusually strong and intra-day variable CP 
has been observed for the BL Lac object PKS 1519—273, with the variability readily explained as 
ISS of a component with up to 3.8% CP, although the origin of the CP itself is unclear (Macquart 
et al. 2000). 

The great sensitivity of SKA will enable the extension of variability modeling to lower lumi¬ 
nosity sources. Monitoring of broad absorption line (BAL) QSOs and other radio quiet objects, 
such as Mrk231, has already revealed variations similar to blazars. Although large scale jets and 
lobes do not develop in these objects, likely due to local environmental conditions—^perhaps the 
presence of the BAL wind, highly energetic flares can still occur close to the core (Reynolds et al. 
2009, 2013). However, with few examples, it is difficult to establish if the flaring is a common 
property in so-called radio quiet AGN. With SKA, the number of low luminosity AGN monitored 
will be extremely large, allowing detailed studies of the variability and comparison to blazars. This 
will give insights into the jet formation processes and bridge the gap between variations in jetted 
systems around stellar-mass black holes and the most extreme blazars. 

Another interesting group of fainter AGN that have recently received a lot of attention are the 
Narrow-Line Seyfert I Galaxies (NLSl). Only 7% of NLSls are formally radio loud (Komossa 
et al. 2006), but interestingly the radio loud ones show many similar properties to blazars, for 
example, several of them are detected in gamma rays by Fermi (Abdo et al. 2009). They also show 
fast variations in their radio light curves, similar to blazars (Foschini et al. 2011; D’Ammando 
et al. 2012). With SKA we will be able to probe the variability in both radio loud and radio quiet 
NLSl objects, shedding light on why some of them are radio loud and others not. If no significant 
variability is detected, it implies that the sources are not Doppler beamed and are viewed at larger 
angles to the line of sight. 

In general, the large number of AGN observed by SKA will allow statistical analysis of AGN 
variability over large populations. By calculating the power spectral densities of the light curves, 
we can gain understanding of the nature of the variations. In cases where significant polarization 
is detected, we obtain information about the underlying magnetic field structure in these objects. 
For example, if the polarization angle is stable over a long time period, it implies that large-scale 
ordered magnetic fields play a role (Aller et al. 2003). The strength of SKA will be the ability to 
study the variations in a large number of sources over a large luminosity range, shedding light on 
AGN unification. 

3.4 Use case: synoptic surveys with outer antennas of SKAl-MID 

It is anticipated that a number of SKAl-MID projects will not make use of the longest available 
baselines. Using the outer antennas to conduct a monitoring programme while the inner antennas 
conduct other observational programmes will be an extremely valuable operational mode. The 
outer antennas could be used for monitoring sufficiently bright targets in order to study their ISS 


15 




Time domain studies of AGN 


Hayley Bignall 


and intrinsic variability, as well as searching for ESEs (and other transient events). Elexible sub¬ 
arraying to allow simultaneous observations at multiple frequencies, and/or of multiple targets, 
would be an important capability. The value of sub-arraying for time domain work was successfully 
demonstrated by the MASIV Survey, which used the VEA split into 5 subarrays to monitor ~ 700 
fields with two-hourly cadence (Eovell et al. 2003). 

3.5 Capabilities of SKAl with 50% reduced sensitivity 

The early science phase of SKA deployment, although not reaching full SKAl sensitivity, 
nevertheless represents a huge leap over existing surveys, and a transition between those possible 
with MWA, MeerKAT, and ASKAP, and the full SKAl. The main effect of 50% reduced sensitivity 
relative to SKAl on the surveys discussed here would be a decrease of a factor ~ 3 in the number 
of sources that could be monitored to a given flux density precision. Eor programmes where sub¬ 
arraying is the desired observing mode, the early science phase would simply use fewer sub-arrays. 
Eimited instantaneous (m,v) coverage of the early roll-out phase may pose problems for the type of 
snap-shot monitoring observations discussed here; in any case an accurate sky model to sufficient 
depth will be required to search for variability. 

3.6 Prospects for SKA2 

Current VEBI surveys reach milliarcsecond resolutions, corresponding to parsec scales in 
AGN. Resolution of the baseline SKA-MID is an order of magnitude or so poorer. Including 
sufficient sensitivity on long baselines as part of SKA2 (see Paragi et al. 2014) would enable the 
expansion of existing monitoring programs to much larger numbers of sources. Polarization ca¬ 
pabilities will be particularly important in enabling identification and tracking of jet components, 
characterizing changes in magnetic fields and their connection to flares, and shedding light on jet 
composition and the interaction of jets with the ambient medium. Simultaneous monitoring of large 
numbers of sources in conjunction with high-energy all-sky instruments will enable us to constrain 
the dissipation of energy during rare AGN outbursts. 

One of the limitations of AGN monitoring at high frequencies is the decrease of survey speed 
with the square of the frequency, assuming single pixel feeds, constant sensitivity, and flat spectrum 
sources. A factor 20 increase in field of view with SKA2-MID would enable survey speeds at 
6 GHz comparable to those with SKAl-MID at 1.4 GHz, enabling the kind of surveys discussed 
above to be undertaken at these higher frequencies, where intrinsic variability becomes increasingly 
important. 

4. Summary 

Variability of AGN provides a powerful probe of the physics of accretion and outflow, as 
well as the dynamics and structure of our own Galaxy’s interstellar medium. The huge leap in 
survey speed engendered by SKA will enable a precision synoptic view of variability for hundreds 
of thousands of AGN, an archive of variability statistics for comparison to ESST and other large 
surveys, and a huge expansion into unexplored parameter space as the study of radio AGN truly 
enters the wide-field time-domain era. These data will allow us to move from the phenomenology 
and statistics of radio variability to understanding the physical mechanisms responsible for the 


16 




Time domain studies of AGN 


Hayley Bignall 


large range of observed variability behaviour. Key eapabilities for this work are broad frequeney 
eoverage—it is partieularly desirable to inelude the 4.6-13.8 GHz Band 5 on SKA-MID; flexible, 
independent sub-arrays whieh ean eover different frequeney bands; and high angular resolution 
eapability. 
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